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The paper considers heat transport in diatomic one-dimensional lattices, containing equal amounts 
of particles with different masses. Ordering of the particles in the chain is governed by single 
correlation parameter - the probability for two neighboring particles to have the same mass. As 
this parameter grows from zero to unity, the structure of the chain varies from regular staggering 
chain to completely random configuration, and then - to very long clusters of particles with equal 
masses. Therefore, this correlation parameter allows a control of typical cluster size in the chain. In 
order to explore different regimes of the heat transport, two interatomic potentials are considered. 
The first one is an infinite potential wall, corresponding to instantaneous elastic collisions between 
the neighboring particles. In homogeneous chains such interaction leads to an anomalous heat 
transport. The other one is classical Lennard-Jones interatomic potential, which leads to a normal 
heat transport. The simulations demonstrate that the correlated disorder of the particle arrangement 
does not change the convergence properties of the heat conduction coefficient, but essentially modifies 
its value. For the collision potential, one observes essential growth of the coefficient for fixed chain 
length as the limit of large homogeneous clusters is approached. The thermal transport in these 
models remains superdiffusive. In the Lennard-Jones chain the effect of correlation appears to be not 
monotonous in the limit of low temperatures. This behavior stems from the competition between 
formation of long clusters mentioned above, and Anderson localization close to the staggering ordered 
state. 


I. INTRODUCTION 

Heat conductivity in one-dimensional (ID) lattices is a 
well-known classical problem related to the microscopic 
foundation of Fourier’s law. The problem started from 
the famous work of Fermi, Pasta, and Ulam (FPU) [1], 
where an abnormal process of heat transfer was initially 
revealed. Numerous aspect of the problem were widely 
addressed over last two decades [2-4]. It was established 
that mere nonlinearity of the interparticle interactions in 
one-dimensional models is insufficient for convergence of 
the heat conduction coefficient in the thermodynamical 
limit. Recently it was stated that if the potential of inter¬ 
action is bounded (allows dissociation of the neighboring 
particles), then the heat conductivity converges [5, 6]. In 
the same time, it is still not clear whether this condition 
is necessary and what is relationship between necessary 
and sufficient requirements for the interatomic potential, 
and dimensionality of the problem. 

In this paper we address other system properties, 
which directly effect the heat transport - a homogene¬ 
ity and an ordering in the lattice. Modifications of the 
heat transport caused by inhomogeneities in the mass or 
potential distribution were explored starting from the pi¬ 
oneering work on harmonic system with different masses 
[7]. Later works considered isotopic disorder [8], har¬ 
monic [9] and annharmonic random chains [10, 11]. How¬ 
ever, main attention of the studies on the nonhomoge- 
neous systems has been attributed to chains with stag¬ 
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gering [12, 13] or randomly distributed masses. In par¬ 
ticular, numerous papers simulated the heat transport in 
one-dimensional diatomic hard-point gas. This system 
is especially interesting, since the case of equal masses 
corresponds to completely integrable linear homogeneous 
billiard (even its version with external on-site potential 
is integrable, see, e.g. [14]). So, the staggering masses 
constitute a perturbation of the integrable case, and the 
question is whether this perturbation will lead to nor¬ 
mal heat transport. Numerous numeric works on the 
diatomic billiard with staggering particles led to a con¬ 
clusion that the heat conduction coefficient k in this sys¬ 
tem diverges in the thermodynamical limit as k ~ N a 
for certain a > 0 [15-18]. More detailed recent explo¬ 
ration of this system [19] demonstrated, that when the 
mass ratio is slightly different from one, it is not possi¬ 
ble to exclude normal heat conduction over longer and 
longer sizes as the integrable limit is approached. In the 
same time, it was conjectured [19] that also in this case 
the heat conductivity diverges, but for longer chains than 
are available for simulation. 

In this paper we also consider the thermal transport 
in the chain models, which include the particles with two 
different masses, but do that in more general setting. Key 
difference from the previous studies is that the ordering of 
these particles in the chains is neither perfectly periodic 
(homogeneous or staggering) nor completely random. In¬ 
stead, the mass distribution in the chain is characterized 
by the correlated disorder. More exactly, the masses of 
neighboring particles are equal with probability p m , and 
0 < Pm < 1. Increase of this parameter favors clustering 
in the chain. All previously studied cases (complete order 
and uncorrelated disorder) correspond to special values of 
p m . The effect of the correlated disorder is studied for the 
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system of rigid particles with abnormal heat conductivity, 
and also for a model with convergent heat conductivity - 
Lennard-Jones (LJ) chain. It is demonstrated that in any 
of these systems the correlated disorder does not modify 
the convergence properties of the heat conduction coeffi¬ 
cient in the thermodynamical limit, but strongly effected 
the quantitative characteristics of the heat transport. 


II. DIATOMIC GAS OF RIGID PARTICLES 
WITH RANDOM MASS DISTRIBUTION 

Let us consider the one-dimensional chain of rigid rods 
with size d > 0, with masses which can take only two 
values Mi = 1 and M 2 = m > 1 (dimensionless param¬ 
eter m is interpreted as the mass ratio). To be specific, 
we choose the length of the rods as d = 0.1 and average 
distance between their centers a = 1. In other terms, 
numeric density of the rods is adopted to be unit. The 
rods are numbered in ascending order of coordinates of 
their centers. Hamiltonian of the system is expressed in 
the following form: 

H = Y,iSr +v(Xn+ (!) 

n n 

Here m n is a mass of the n-th rod, x n - coordinate of its 
center, p n = m n x n . The interaction of absolutely rigid 
particles is described by the following hard-core poten¬ 
tial: 

V(r) = 00 if r < d, V(r) = 0 if r > d. (2) 

Potential function (2) describes instantaneous elastic col¬ 
lisions of the rods. In this one-dimensional chain, only 
neighboring rods can collide, and collision of the rod num¬ 
ber n with the rod number n + 1 occurs if the distance 
between their centers x n +± — x n = d. As it is well-known, 
if before the collision velocities of the rods are v n = x n 
and v n +i = x n +i, then after the collision the velocities 
will take the following values: 

v' n = [2m n+1 v n+1 + (m„ - m n+1 )v n ]/(m n + m n+1 ) 
v n+i = [2m n i>„ + (m n+ 1 - m n )v n+1 }/ (m n + m n+ 1 ). 

Between the collisions, the rods move as free particles. 

In considered model, each rod can have the mass 
m n = 1 or the mass m n = m with equal concentration. 
To describe the correlation between the masses of neigh¬ 
boring rods, we define additional parameter 0 < p m < 1, 
which denotes a probability for the neighboring rod to 
have the same mass. The case p m = 0 corresponds to the 
chain with alternating masses, the case p m = 0.5 corre¬ 
sponds to completely random distribution of the masses 
in the chain (lack of correlation between the neighbors). 
The higher value of p m , the longer clusters of particles 
with equal equal masses are expected in the chain. Av¬ 
erage length of such homogeneous clusters is estimated 


as N p ~ 1/(1 — p m ). In the limit p m —> 1 the chain be¬ 
comes almost homogeneous, in other terms, it contains 
the homogeneous clusters of diverging length N p —)• 00 ). 

For simulation of the heat transport in this model with 
N rods, we include the interaction of terminal rods with 
boundary thermostats. The rod with n = 1 interacts 
with thermostat of temperature T+, the rod number N 
- with thermostat of temperature T_. Interaction of the 
first rod with the thermostat occurs when X\ = d/2. In 
this moment, the velocity of this rod is re-ascribed to 
V\ > 0; the latter is random with Maxwell distribution 

P(v ) = (\v\m/T) exp(— v 2 m/2T) 

with mass m = m\ and temperature T = T+. Simi¬ 
larly, the rod N interacts with Maxwell thermostat when 
xn = N — d/2. In this moment, the velocity of this is 
re-ascribed to random value vn < 0 with Maxwell distri¬ 
bution with mass m = tun and temperature T = T_. 

As it was mentioned above, the rods interact with 
the thermostats only when they collide with the bound¬ 
aries. In the moment of collision t = tj the thermo¬ 
stat changes the energy of the terminal rod by value 
A Ei(tj) = rrii[vf(tj + 0) — vf(tj — 0)]/2, where i = l,N. 
If over time interval [0, t] there were N t collisions of 
the terminal rod with the boundary in time instances 
1 e [0,£]), then the average work done by the 
thermostat is expressed as 

1 Nt 

3i(t) = 

j =1 

and its average power Ji = Hindoo j,; {/). 

For simulation of the heat transport in the system we 
choose the following initial conditions: 

*n(0) = n - 1/2, x n (0)=t; n , n=l,2, 

Here v n is a random value with Maxwell distribution 
P(v) = y/m n /2TrT exp[— m n v 2 /2T], where T = (T + + 

T-)/ 2- 

Following paper [19], the temperature of the left 
boundary is set to T + = 6, and of the right boundary - to 
T_ = 4. Then, long time dynamics of the system is sim¬ 
ulated. It should be mentioned that this dynamics does 
not depend on the absolute values of the temperatures, 
but only on the ratio T + /T_. After initial transient and 
formation of stationary heat flux, average powers of the 
thermostats Ji, Jn are computed. The heat flux in the 
system should satisfy relationship J = J\ = — Jn, which 
can serve also as a test for validity of the numeric pro¬ 
cedure. In the simulations these equations were obeyed 
with very high accuracy. Distribution of the local tem¬ 
perature in the chain is determined as T n = (i^(£)) t . 

Simulation of the heat transport demonstrates that the 
collisions of the rods with different mass bring about for¬ 
mation of a linear thermal gradient in the chain, see Fig¬ 
ure 1. The effect of thermal resistance reveals itself at 
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FIG. 1: (Color online) Distribution of temperature T n along 
the chain of length N = 160. The mass ratio is m = 1.4, 
parameter of correlation p m = 0. The dashed line represents 
the linear temperature gradient, used for computation of the 
heat conduction coefficient k(N). 



the ends of the chain. As a result of this effect it turns 
out that the temperature of the leftmost rod T\ < T + , 
and of the rightmost rod Tn > T_. This discrepancy 
decreases with an increase of the chain length. To avoid 
ambiguities related to this thermal resistance, value of 
the heat conduction coefficient was calculated in accor¬ 
dance with actual temperatures of the terminal rods with 
the following formula: 

k(N) = JN/(T 1 -T n ). (3) 

Dependence of the heat conduction coefficient re on the 
chain length N for different values of the mass ratio m 
and mass correlation parameter p m is presented in Fig¬ 
ure 2. One can learn from this Figure, that for the case 
Pm = 0 (perfectly ordered diatomic chain) the depen¬ 
dence re(JV) is similar to results presented previously in 
paper [19]. For large value of the mass ratio to = 3 
the heat conduction coefficient increases approximately 
as TV 0 - 3 . However, for small mass ratio m = 1.1 for 
N ~ 10 4 the growth of the heat conduction coefficient 
slows down. It is extremely difficult to check the behav¬ 
ior of the coefficient for essentially longer chains. 

Increase of the mass correlation parameter p m for fixed 
mass ratio and chain length leads to an increase of re, as 
it is demonstrated in Figures 2, 3. In the limit p m —> 1 
one obtains re(TV) —> oo. The reason is that the modifi¬ 
cation of momenta of the individual rods is possible if the 
colliding particles have different masses. Average length 
of clusters of the rods with the same mass N p grows to¬ 
gether with the mass correlation parameter p m , and rela¬ 
tive number of the scattering events collisions of the rods 
with different mass decreases. So, one can expect that 
for large N the heat conduction coefficient re will be pro¬ 
portional to the average length of the cluster with equal 
masses, which can be estimated as N p ~ 1/(1 —p m ). As 
it follows from Figure 3, in the case of large mass ratio 
to = 3 we observe the divergence of the heat conduction 


FIG. 2: (Color online) Dependence of the heat conduction 
coefficient re on the length of chain N for mass ratios m = 1.1, 
1.4, 3 and for three values of the mass correlation parameter: 
Pm = 0 (blue curves 1, 2, 3), p m = 0.5 (red curves 4, 5, 6), 
Pm = 0.75 (green curves 7, 8, 9). 



FIG. 3: (Color online) Dependence of the heat conduction 
coefficient re on the mass correlation parameter p m for chain 
length N = 1280 (curves 1, 2, 3), N = 2560 (curves 4, 5,6) 
and mass ratios m = 1.1, 1.4, 3. Straight line corresponds to 
expression re = 3/(1 — p m )- 


coefficient re ~ 1/(1 —p m ). For smaller values of the mass 
ratio, no clear power law is observed. For m = 1.4 it is 
still possible to argue approaching the power law, but in 
the case of to = 1.1 the crossover is much slower, and 
simulation of longer chains is required. 

The convergence behavior of re with increase of p m 
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FIG. 4: Transport of initial energy pulse in the chain with 
N = 20000 particles after initial thermal excitation at the left 
end. (a) m = 1, time step dt = 144; (b) m = 1.4, dt = 3600; 
(c) m = 3, dt — 3960, p m = 0. Energy distribution in the 
chain E(x) is plotted for multiples of dt. 


strongly resembles the situation with m —> 1. 

As it was mentioned before, our results for p m = 0 
reproduce simulations presented in [19]. There the au¬ 
thors have mentioned that on the basis of the simulation 
results it is not possible to reject possibility of conver¬ 
gence of the heat conduction coefficient for the case of 
small mass ratios. From Figure 2 one can observe that 
for larger values of the mass correlation parameter the 
situation remains qualitatively the same. The authors 
of [19] also conjectured that the observed ’’convergence¬ 
like” behavior is transient, and the heat conduction co¬ 
efficient should attain the divergence regime for higher 
values of N. Results presented in Figure 3 give a clue 
that the divergence of k in the limit of p m —> 1 exhibits 
similar very long crossovers. In the next Section, we at¬ 
tempt to achieve numeric evidence on divergence of the 
heat conduction coefficient in the ’’crossover” situations 
by simulating the dynamics of temperature spreading in 
the chain. 


III. DIFFUSION OF ENERGY PULSES IN THE 
CHAIN WITH RANDOM MASS DISTRIBUTION 

For simulation of the energy diffusion the chain com¬ 
prising N = 20000 particles (rods) is considered. The 
thermostats are removed and the terminal rods elasti¬ 
cally interact with the walls. Initially all the rods are 
randomly placed in the interval (d, N — d) and 60 rods 
at the left end of the system are thermalized with aver¬ 
age temperature T in . The rest of the rods (for n > 60) 
are thermalized with relatively small average tempera¬ 
ture To < Ti n ■ Then, the diffusion of energy along the 
chain is simulated. To improve the accuracy, all results 
of such simulations are averaged over 10 4 independent 
realizations of the initial temperature distribution. For 
every time instance one obtains energy distribution in 
the chain 

E n (t) = ^2 mi±f(t)/2. 

i:xi£(n—l,n\ 

Let us at first consider energy diffusion in initially non- 
thermalized chain. For this matter we adopt the value 
Ti n = 10 as the initial temperature of the left end of the 
chain, and the temperature of the rest part of the chain 
is T 0 = 0. Results of these simulations are presented in 
Figure 4. The character of the thermal diffusion depends 
on the value of mass ratio m. For m = 1 one observes 
completely ballistic flow all energy leaves the left end of 
the chain, and the heat pulse becomes wider and lower, 
since the initial pulse contains particles with different 
velocities. For m > 1 the flow includes both ballistic and 
diffusive components, and some residual energy remains 
at the left end of the system. 

It is common to describe evolution of the temperature 
distribution in the chain with the help of the following 
’’mean quadratic” deflection from the initial state (con¬ 
centration at the left end of the chain): 

N 

a 2 = J^(n - l/2) 2 e n . (4) 

n —1 

Here {e n = (E n — T 0 )/2T}() r =1 characterizes the energy 
distribution in given time instant, E = J22 = i{E n — -Eo) is 
a total constant energy of the chain, Eo = T 0 /2 - energy 
level of initially ” cold” part of the chain. 

Numeric simulation demonstrates that for non- 
thermalized chain (To = 0) the mean quadratic deflection 
at large times asymptotically obeys power law: (a 2 ) ~ 
when t /* oo. These results are illustrated in Figure 5. 

The works [20-23] argue that in systems with collision 
dynamics, where the particles undergo Levy flights, the 
exponent /3, which describes the mean square displace¬ 
ment of particles (x 2 (t)) ~ , should be related with a 

through equation 

a = /3 — 1. (5) 

This equation implies that normal diffusion (/3 = 1) leads 
to the normal (non-divergent) heat conductivity (a = 0). 
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FIG. 5: (Color online) Time dependence of the measure of 
energy spreading over the chain a 2 (double logarithmic scale) 
for the set of parameters: m = 1 (curve 1); m = 1.05 and 
Pm = 0 (curve 2); m = 1.1 and p m = 0, 0.5, 0.75 (curves 
3, 4, 5 respectively); p m = 0 and m = 2, 3 (curves 6, 7 
respectively). Dashed lines correspond to power laws for 
the values of exponents /3 = 2, 1 and 1.31 (curves 8, 9 and 10, 
11 respectively). The temperature of the chain is To = 0. 

In the same time, superdiffusion (j3 > 1) corresponds to 
divergent (a > 0) heat conductivity. 

Numeric results in Figure 5 for homogeneous chain 
with m = 1 yield value /3 = 2, which corresponds to 
purely ballistic energy transport (a = 1). For mass ra¬ 
tios m = 1.05, 1.1, 2 and 3 we obtain f3 = 1.31, which is 
related to the superdiffusion of energy in non-thermalized 
diatomic chain with mass ratio m > 1. The correlation 
in mass distribution along the chain does not lead to a 
change in the type of diffusion. So, for m = 1.1 the value 
of the exponent /3 does not depend on the value p m < 1. 
For m = 3 the measured exponent of divergence of the 
heat conductivity is a = 0.3, such that for all the values 
p m there is a good agreement with the relation (5). 

Let us consider now diffusion of energy in a thermal- 
ized chain with To > 0. Figure 6 illustrates that dy¬ 
namics of the energy pulse significantly depends on the 
temperature of the chain. An increase of the temperature 
Tq leads to an increase in the rate of propagation of the 
heat pulse along the chain. The exponent of the energy 
diffusion /3 increases with an increase in the mean tem¬ 
perature To - see Fig. 7. So, in the diatomic chain with 
staggering masses ( p m = 0) for mass ratio m = 2 with 
To = 0 we obtain the energy diffusion exponent /3 = 1.31, 
for T 0 = 0.1 P = 1.69, for T 0 = 0.2, /3 = 1.73, and for 
T 0 = 1.77 the value of the exponent /3 = 1.77. We are 
not aware of any significant dependence of a on average 
temperature (or, in this collisional model, on ratio of the 




FIG. 6: Transport of initial energy pulse in the chain with 
N = 20000 particles after initial thermal excitation at the left 
end Tin = 10 + To for temperature of the chain (a) To = 0, 
time step dt = 3600; (b) To = 0.2, dt = 1500 and (c) To = 0.4, 
dt = 1320. Mass correlation parameter is p m = 0, mass ratio 
m = 2. Energy distribution in the chain E(x) is plotted for 
multiples of dt. 


temperatures on the cold and the hot ends of the chain). 
Thus, the obtained results suggest that the relationship 
between a and /3 (5) is satisfied only for To = 0. Both 
in the cases To > 0 and To = 0 the value of the diffu¬ 
sion exponent /? does not exhibit any visible dependence 
on the correlation of the mass distribution of the rods in 
the chain. The correlation parameter p m only effects the 
time at which the system achieves the regime of steady 
superdiffusion - see Fig. 8. 

The results of the numerical modeling in a random 
diatomic chain imply that for mass ratio m > 1 and for 
any correlation of the mass distribution of the chain, i.e. 
for any value of the parameter 0 < p m < 1, the heat 
transport in the diatomic chain is superdiffusive for all 
studied values of mass ratio and correlation parameter. 
These results include also the cases, in which it was not 
possible to get clear evidence of divergence of the heat 
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FIG. 7: (Color online) Time dependence of the measure of 
energy spreading over the chain <r 2 (double logarithmic scale) 
for temperature of the chain To = 0, 0.1, 0.2 and 0.4 (curve 
1, 2, 3 and 4). Mass ratio m = 2, mass correlation parameter 
p m = 0. Dashed straight lines correspond to power laws for 
the values of exponents /3 = 1.31, 1.69, 1.73 and 1.77 (curves 
1, 2, 3 and 4 respectively). 


FIG. 8: (Color online) Time dependence of the measure of 
energy spreading over the chain a 2 (double logarithmic scale) 
for the chain with mass correlation parameter p m = 0, 0.5 
and 0.75 (curve 1, 2 and 3). Mass ratio m = 1.1, temperature 
of the chain To = 0.2. Dashed straight line corresponds to 
power law t 1 ' 75 . 


conduction coefficient in the equilibrium simulations. 


IV. HEAT TRANSPORT IN RANDOM 
DIATOMIC LENNARD-JONES CHAIN 


To explore the effect of mass clustering in the diatomic 
lattice with convergent heat conductivity, we consider a 
chain with common Lennard-Jones (LJ) interatomic po¬ 
tential: 


V(r) = 4e[(er/r) 6 — 1/2] 2 . (6) 

The following values of parameters are used: e = 1/72, 
a = 2 -1 / 6 . For this choice of the parameters, the po¬ 
tential has a minimum for the distance between particles 
ro = 1 and stiffness in the minimum k = V"(ro) = 1. 
Consequently, the soun d ve locity in the homogeneous 
(to : 1) chain vq = ro^/k/m = 1. 

Recently it was established that the homogeneous 
Lennard-Jones chain has convergent heat conductivity 
[5, 6]. Efficient phonon scattering exists due to thermally 
activated large elongations of the interparticle bonds; 
possibility of such large elongations stems from bound¬ 
edness of the LJ potential [6]. All simulations described 
below confirm the convergence of the heat conduction in 
this system, although this issue is still debated [24] . One 
might expect that inhomogeneities in the chain enhance 
the phonon scattering and, as a consequence, decrease 


the heat conductivity. Therefore for to > 1 the heat con¬ 
ductivity of the LJ chain should be convergent, as in the 
homogeneous case. 

To check this idea, the heat transport in random 
LJ chain with end regions attached to Langevin ther¬ 
mostats is simulated. The method of simulation is de¬ 
scribed in details in paper [5]. The LJ chain contains 
N = 80 + 10 x 2 l , l = 1,2,..., 10 particles with fixed 
boundary conditions. Thee leftmost 40 particles are at¬ 
tached to the Langevin thermostat with temperature T + 
and 40 rightmost particles - to the Langevin thermostat 
with temperature T_, where T± = (1 ± 0.1)T. Here T 
is the average temperature of the chain. Local tempera¬ 
ture distribution in the chain is defined as T n = ( x 2 ), the 
heat flux - as J = (J n ), where J n = —x n V'(x n +i—x n ) is 
a local instantaneous flux of energy through site number 
n. In internal fragment of the chain 40 < n < N — 40 
one observes constant heat flux and approximately linear 
thermal profile. Thus, it is possible to evaluate the heat 
conduction coefficient as /t(iVj) = JNi/(Tn — Tn- 40 ), 
Ni = N — 81 is the length of the chain fragment between 
the thermostats. 

Thermal conductivity can also be found with the help 
of Green-Kubo relation [25] 


k = 


lim 

t—too 


1 


lim 

AT^oo NT 2 



( 7 ) 


where c{t ) = ( J s (t)J s (t — t)) T autocorrelation function 
of the heat flux; the latter is defined as J s {t) = J2 n J n (t). 
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FIG. 9: (Color online) Dependence of the thermal conduc¬ 
tivity k on the length N of the diatomic LJ chain for tem¬ 
perature (a) T = 0.0005 and (b) T = 0.2. The mass ratio 
m = 1.1. The curves represent the results for three values of 
the mass correlation parameter p m = 0, 0.5 and 0.75 (curves 
1, 2 and 3 respectively). Horizontal lines correspond to the 
values of the heat conduction coefficient obtained with the 
help of Green-Kubo relation. 


To find the autocorrelation function c(f), we simulate 
a cyclic chain with N = 10 4 units, initially completely 
immersed in the Langevin thermostat with temperature 
T. After thermal equilibration with the thermostat, the 
latter is removed and Hamiltonian dynamics of the chain 
is simulated. To increase the accuracy of the autocorre¬ 
lation function measurement, its value is averaged over 
10 4 independent realizations of the initial chain thermal- 
ization. 

In the case of small temperatures T <C e the heat 
transport occurs due to propagation of weakly interact¬ 
ing linear waves (as it is commonly accepted, we use the 
term ’’phonons” despite purely classical character of our 
model). If the temperatures arc large, T e only hard 
cores of the interatomic potentials have significant effect 


on the heat transport. So, the LJ chain can be consid¬ 
ered as a set of almost completely rigid particles, and the 
heat transport occurs due to the collisions between these 
particles. Since the heat transport mechanisms for high 
and low temperatures are rather different, it is natural 
to expect different effect of the random mass distribution 
in these two cases. 

First, let us consider the heat transport in the LJ 
chain for relatively low temperature T = 0.0005. In this 
limit the system becomes very close to the chain with 
parabolic interatomic potential. Therefore, the diatomic 
LJ chain with staggering masses (mass correlation pa¬ 
rameter p m = 0) becomes close to well-known homoge¬ 
neous diatomic chain with its typical acoustical and opti¬ 
cal phonon branches. In the same time, for random chain 
with mass correlation parameter 0 < p m < 1 it is well- 
known that some vibration modes become localized (An¬ 
derson localization). So, one can expect that as p m grows 
above zero and the homogeneous diatomic chain becomes 
randomized, the heat conductivity of the chain should 
decrease. In the same time, as p m approaches unity, the 
chain becomes closer to the monoatonric case, and the 
heat conductivity should increase. Therefore, contrary to 
the chain of rigid particles explored above, one can quali¬ 
tatively predict non-monotonous dependence of the heat 
conduction coefficient on the mass correlation parameter 
p m for fixed temperature and mass ration. Moreover, it 
is anticipated that this effect will be more pronounced 
for higher mass ratio. 

Dependence of the heat conduction coefficient k on the 
chain length N for different values of p m is presented in 
Fig. 9. For the mass ratio close to unity m = 1.1, the 
increase of p m leads to relatively mild decrease of the 
heat conductivity only for short chains with N < 10 3 , as 
it is illustrated in Figure 9 (a). Possible explanation is 
that the chain is rather close to the homogeneity, thus 
the Anderson localization is less pronounced and other 
effects (such as grow of expected length of the homoge¬ 
neous clusters) can compete and blur the effect of local¬ 
ization. For longer chains, the heat transport is governed 
by phonons with large wavelength, which are less sensi¬ 
tive to small inhomogeneities of the mass distribution. 

For high temperature T = 0.2 the chain dynamics is 
primarily governed by collisions between hard repulsive 
cores of the particles. It is possible to say that in this 
case the LJ chain is close ID diatomic gas model, but 
not completely equivalent to it. The crucial difference 
is that for any finite values of the temperature the time 
of collisions is finite. Therefore, contrary to the ID gas 
considered above, triple collisions of the particles have 
a nonzero probability. In recent paper [6] it is demon¬ 
strated that such triple collisions provide efficient scat¬ 
tering mechanism for energy transport in the system and 
promote convergence of the heat conduction coefficient. 
Consequently, one should expect a monotonous growth of 
the heat conductivity with increasing p m similarly to the 
diatomic gas. Numeric simulation presented in Figure 9 











FIG. 10: (Color online) Dependence of the thermal conduc¬ 
tivity k on mass correlation coefficient p m in the LJ chain 
with mass ratio m = 1.1 and temperatures T = 0.0005 and 
T = 0.2 (curves 1 and 2) and for m = 2, T = 0.0005 and 
T = 0.01 (curves 3 and 4). 

(b) confirms this similarity (cf. Figure 3). Important 
difference is that in the LJ chain the heat conductivity 
for p m — > 1 remains convergent in the thermodynamical 
limit. 

Due to this convergence, it is reasonable to consider 
a dependence of the limit value of the heat conduc¬ 
tion coefficient on the mass correlation parameter p m . 
This dependence has been analyzed, for the chain with 
N = 20000 particles for different temperatures and mass 
ratios with the help of Green-Kubo relation (7). Method¬ 
ology of this simulation has been described above, and 
the results are presented in Figure 10. 

From Figure 10 one can learn that the character of 
dependence k{ p m ) substantially depends on the temper¬ 
ature and mass ratio. Common feature of all these de¬ 
pendencies is a rapid growth of the heat conductivity as 
p m —>■ 1. For relatively low temperature T = 0.0005 
and for both considered mass ration m = 1.1 and m = 2 
these dependencies are not monotonous. In particular, 
for m = 1.1 the heat conductivity achieves a minimum 
for p m = 0.2 and for m = 2 the minimum is achieved for 
p m = 0.75 (see Fig. 10, curves 1 and 3). As it was men¬ 
tioned above, this non-monotonicity may be attributed 
to Anderson localization of oscillatory states in the chain 
with random structure - two terminal states with p m = 0 
and p m = 1 are ordered. This effect is much more pro¬ 
nounced for higher mass ratio m = 2. 

For the case of relatively high temperatures this phe¬ 
nomenon of non-monotonicity disappears. Namely, for 
mass ratio m = 2 and temperature T = 0.01 the 


heat conductivity is almost constant over the interval 
0 < Pm < 0.9 and rapidly grows for p m > 0.9 (see Fig¬ 
ure 10, curve 4). For mass ratio m = 1.1 and temperature 
T = 0.2 the heat conductivity exhibits a monotonous 
growth over the whole interval 0 < p m < 1 (see Fig¬ 
ure 10, curve 2). This behavior is another illustration 
of primarily collisional dynamics of the LJ chain in the 
limit of high temperatures - it behaves similarly to the 
gas of rigid particles, cf. Figure 3. It is important to 
observe that in this regime there is no significant phys¬ 
ical difference between the ordered staggered state with 
Pm = 0 and disordered states with small values of p m . 
The reason is that the collision mechanism is local and 
not sensitive to the short-range order or disorder. 


V. CONCLUDING REMARKS 

The first conclusion to mention here is that the inter¬ 
nal correlations in the particle ordering do not modify 
the convergence properties of the heat conduction co¬ 
efficient. In the model of colliding rods, superdiffusive 
transport has been detected for all explored values of the 
mass ratio an the correlation coefficient, including those 
for which direct equilibrium simulation of the heat trans¬ 
port is beyond current computational capabilities. It is 
worth while mentioning, that the superdiffusion expo¬ 
nent turned out to be dependent on initial heating of the 
chain beyond the initial thermal pulse. 

In the same time, the correlations in the mass disorder 
have significant effect on the value of the heat conduc¬ 
tivity. This effect reveals itself through two competing 
factors - growth of the clusters of particles with the same 
mass, and frustration of the perfect disordered state. In 
the case of purely collisional dynamics, the latter factor 
plays no role, and the heat conductivity of the chain with 
given length monotonously increases with increase of the 
correlation parameter. In the case of relatively small tem¬ 
peratures in LJ chain both factors are significant, and the 
dependence is not monotonous due to Anderson localiza¬ 
tion as the system departs from the ordered staggering 
state. For large temperatures, the quantitative depen¬ 
dence of the heat conduction coefficient on the correla¬ 
tion parameter is similar to that in the model of colliding 
rods; however, the thermal transport appears to remain 
convergent. 
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